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Abstract 

This is the first of a series of papers in which we examine the con- 
straints of sphericaUy symmetric general relativity with one asymptot- 
ically flat region. Our approach is manifestly invariant under spatial 
diffeomorphisms, exploiting both traditional metric variables as well as 
the optical scalar variables introduced recently in this context. With 
respect to the latter variables, there exist two linear combinations of 
the Hamiltonian and momentum constraints one of which is obtained 
from the other by time reversal. Boundary conditions on the spher- 
ically symmetric three-geometries and extrinsic curvature tensors are 
discussed. We introduce a one-parameter family of foliations of space- 
time involving a linear combination of the two scalars characterizing a 
spherically symmetric extrinsic curvature tensor. We can exploit this 
gauge to express one of these scalars in terms of the other and thereby 
solve the radial momentum constraint uniquely in terms of the radial 
current. The values of the parameter yielding potentially globally reg- 
ular gauges correspond to the vanishing of a timelike vector in the 
superspace of spherically symmetric geometries. We deflne a quasi- 
local mass (QLM) on spheres of flxed proper radius which provides 
observables of the theory. When the constraints are satisfled the QLM 
can be expressed as a volume integral over the sources and is positive. 
We provide two proofs of the positivity of the QLM. If the dominant 
energy condition (DEC) and the constraints are satisfled positivity can 
be established in a manifestly gauge invariant way. This is most easily 
achieved exploiting the optical scalars. In the second proof we specify 
the foliation. The payoff is that the weak energy condition replaces 
the DEC and the Hamiltonian constraint replaces the full constraints. 
Underpinning this proof is a bound on the derivative of the circum- 
ferential radius of the geometry with respect to its proper radius. We 
show that, when the DEC is satisfled, analogous bounds exist on the 
optical scalar variables and, following on from this, on the extrinsic 
curvature tensor. We compare the difference between the values of the 
QLM and the corresponding material energy to prove that a reason- 



able definition of tfie gravitational binding energy is always negative. 
Finally, we summarize our understanding of the constraints in a tenta- 
tive characterization of the configuration space of the theory in terms 
of closed bounded trajectories on the parameter space of the optical 
scalar s. 



I INTRODUCTION 

To identify the independent dynamical degrees of freedom of the gravitational field in 
general relativity it is useful to cast the theory in Hamiltonian form [1] . This means that the 
gravitational field must be viewed, not as a fixed four-dimensional object, but rather as a 
sequence in 'time' of Riemannian three-geometries. Thus we think of specifying some initial 
configuration of sources and gravitational field and letting it evolve. As is well known, this 
initial data cannot be specified arbitrarily, it must satisfy the Einstein constraint equations. 
These constraints only contain the source energy density and momentum density. They do 
not depend on the equations of state. Of course, if we wanted to track the evolution of the 
system we would need to provide a more detailed specification of the sources, including 
these equations of state. 

In this paper we will focus on the solution of the classical constraints and the identifi- 
cation of those features of the theory which depend only on the initial data in the simplified 
setting of spherically symmetric general relativity with our sight set on the quantum theory. 

A remarkable consequence of the diffeomorphism invariance of general relativity is 
that, in a sense, the constraint equations are all there is to the theory. For if the constraints 
are satisfied at all times and the sources are completely specified, then the evolution 
equations follow [2]. For this reason, the solution of the constraints should be viewed 
as much more than a prerequisite to the solution of the dynamical problem. Once a point 
in the classical configuration space {i.e. a solution of the constraints) has been identified, 
its subsequent evolution is implicitly defined. The structure of this space is, of course, 
highly non-trivial. 

The most developed classical approach to the solution of the constraints has been the 
conformal geometry approach pioneered by York and co-workers in the seventies. This was 
very successful in settling formal questions such as the existence and uniqueness of solu- 
tions [3]. However, beyond this formal level, it is extremely difficult to piece together the 
structure of the configuration space of the full theory outside the domain of perturbation 
theory. Unfortunately, this is the framework on which the canonical quantization of the 
theory is based. Thus, until this is done any claims we make about the quantum theory 
must necessarily be taken with a grain of salt. 

One regime in which the problem simplifies, without sacrificing all local dynamical 
degrees of freedom (such as we do in homogeneous relativistic cosmologies), is when the 
geometry as well as the material sources are spherically symmetric [4,5,6,7]. In such a 



system, all the true local dynamical degrees of freedom reside in the sources. There are 
no independent local gravitational dynamical degrees of freedom. The sources, however, 
generate a 'gravitational potential', a kinematical object, which in turn interacts on them. 
The dynamics of matter associated with this potential can be extremely non-trivial, a point 
convincingly demonstrated by the recent controversy created by Choptuik's numerical 
simulations of the collapse of a massless scalar field [8] . 

There are only four topologies compatible with spherically symmetric initial data that 
is defined on a three- manifold. The manifold can be R^, with a regular center and one 
end, just like ordinary flat space; it can be S'^ x R^, with two ends and no center as with 
the spatial slice through extended Schwarzschild spacetime; it can be S^ x S^ which is the 
spherically symmetric torus or it can be S^, the three-sphere [9]. We limit the discussion 
to the first case, i.e., to geometries possessing one asymptotically flat region deferring the 
examination of spherically symmetric inhomogeneous cosmologies and the double ended 
case to future publications. 

The boundary conditions associated with the given topology play an important role. In 
the case we will study, the only boundary condition we need to implement is the regularity 
(or the degree of singularity) at the base of the spatial geometry. Technically this is 
because the Hamiltonian constraint is a singular ordinary differential equation at this point. 
On one hand, this imposes an extraordinary rigidity on the solution, making it unique. 
On the other it provides the mechanism, when the energy density is appropriately large, 
which allows singularities to occur in the geometry. If the material sources are suitably 
localized (as we will always assume) the constraints will automatically steer the geometry 
to asymptotic flatness if no singularity intervenes. In closed cosmologies, the non-singular 
closure of the spatial geometry imposes integrability conditions on the sources it contains. 

The initial data for the gravitational fleld consists of two parts, the intrinsic geometry 
of the three-manifold and an extrinsic curvature tensor which describes how this three- 
manifold is embedded into a four dimensional spacetime. The solution of the constraints 
involves the implementation of gauge conditions. One of these conditions involves the 
speciflcation of how these three-manifolds foliate spacetime. There are two ways of doing 
this within the canonical context; intrinsically, where the foliation is determined by placing 
some restriction on the three- metric, for example, that it be flat; and extrinsically, where 
some condition is placed on the extrinsic curvature, for example, that its trace vanishes 
(the maximal slicing condition). For any choice one must show that it is compatible with 



the constraints and that it can be used as an evolution condition. It has been found that 
doing this extrinsicaUy is invariably better than doing it intrinsically. In view of this we 
will only consider extrinsic slicing conditions. 

The remaining gauge condition concerns the specification of the spatial coordinate 
system. The point of view we will adopt in this paper is that it is not necessary, at least at 
the level of the constraints, to make an explicit spatial coordinate choice. The justification 
for this is the fact that there are two invariant linear measures of the spherically symmetric 
geometry, the circumferential radius R, and the proper radial length i, and the constraints 
come ready cast in terms of derivatives of R with respect to i. 

It is natural that the gauge which fixes the foliation is the gauge which should be 
tackled first. Fix the foliation, then fix coordinates on the hypersurfaces picked out by this 
foliation. Having said this, it is only fair to also point out that the choice of gauge which 
simplifies the solution of the constraints most dramatically is implemented most efficiently 
by inverting this order, exploiting the circumferential radius as the radial coordinate and 
then foliating spacetime by the so called polar gauge* In this gauge, not only does the 
extrinsic curvature quadratic miraculously fall out of the Hamiltonian constraint so that 
it mimics its form at a moment of time symmetry, but the constraint is then also exactly 
solvable. Furthermore, the momentum constraint reduces to an algebraic equation which 
permits the non-vanishing extrinsic curvature component to be determined locally in terms 
of the material current. What is unfortunate is that both the foliation and the spatial 
coordinate system break down catastrophically when the geometry possesses an apparent 
horizon. This corresponds to the vanishing of one or the other of Q±, the divergence of 
the future and past pointing outward directed null rays on a metric two-sphere at fixed 
proper radius [10]. 

For the purpose of examining observable effects in the classical theory it is sufficient to 
truncate the geometry at the horizon if it possesses one, and place appropriate boundary 
conditions there. Even if the formation of the horizon is a consequence of physical processes 
occuring in its interior once formed the details of the interior physics can have no observable 
consequences in the exterior. In the quantum theory, however, we know that we are not 
always at liberty to truncate the theory in this way. 

On one hand, a process like the Hawking effect can be understood in terms of the 
polarization of the vacuum in the exterior neighborhood of the event horizon [11]. In the 



* This is the gauge exploited in refs. [4,5] and [6]. 
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approximation in which the back-reaction on the geometry can be ignored the techniques of 
quantum field theory on a given curved background spacetime apply. What is beyond the 
scope of any approximation which truncates the geometry at the horizon is the prediction 
of the final state of the black hole. 

There are other processes, however, still below the Planck regime, such as tunnelling 
from a configuration without an apparent horizon to a configuration with an apparent 
horizon, in which the existence of classically inaccessible regions of the spatial geometry 
can have dramatic consequences in the quantum theory [12,13]. 

The above processes in quantum gravity are both semi-classical in nature. In the 
Planck regime, however, we do not even possess an unambiguous classical lump to start 
with. Furthermore, the very definition of an apparent horizon involves both the intrin- 
sic and the extrinsic geometry (or equivalently the momentum conjugate to the intrinsic 
geometry) which are not simultaneously observable in the quantum theory. 

If we adhere to a configuration space consisting of metric variables, for the canonical 
quantization of the model we will need to catalog all possible solutions satisfying the 
constraints with or without apparent horizons. The only way to mend the situation to 
accomodate the polar gauge would be to introduce the gauge patch by patch between 
successive horizons. We do not examine this possibility here because it would be almost 
impossible to implement in the quantum theory. 

A canonical change of variables, from the traditional metric phase space variables, 
to the optical tensor variables defined on a foliating sequence of closed two-dimensional 
hypersurfaces embedded in the three-geometry provides an extremely useful alternative 
description of the initial data when the geometry is spherically symmetric. In this case, 
when the two dimensional hypersurfaces are also spherically symmetric, the optical tensors 
reduce to the two scalar quantities, G±. The vanishing of G_|_ corresponds to a future 
apparent horizon, and the vanishing of G_ to a past apparent horizon. Thus, by adopting 
this variable to characterize the configuration space, we sidestep the difficulty inherent in 
the metric variable description of apparent horizons in the quantum theory. These variables 
are a linear combination of intrinsic and extrinsic quantities [14]. Most importantly, is that, 
when cast with respect to the optical scalars, we can replace the Hamiltonian constraint 
and the momentum constraint by a pair of quasi-linear first-order equations, one of which 
is the time reversal of the other [15] and which are entirely equivalent to the original 
constraints. 



In Sect. 3 we return to the metric variables in a search for a globaUy vahd fohation. 
We introduce a one-parameter family of foliations corresponding to the vanishing of some 
linear combination of the two independent scalars characterizing the extrinsic curvature in 
a spherically symmetric geometry. Each such gauge corresponds to a ray in superspace [16]. 
The physically acceptable foliations correspond to timelike directions. Maximal slicing is 
one of these. With the optical variables, this is the natural choice of gauge. However, we 
find that there are other unexpected parameter values possessing attractive features. One 
of the lightlike directions in superspace bounding the valid gauges corresponds to polar 
gauge. In the gauge defined by the other lightlike direction, the Hamiltonian constraint 
also mimicks its form at a moment of time symmetry. As such, it is worth considering more 
closely. Minimal surfaces in this gauge, however, do not coincide with apparent horizons. 
What is more serious, the foliation is not suitably aymptotically fiat. 

One of the most remarkable results in general relativity is the positivity of the Arnowitt-| 
Deser-Misner (ADM) mass, the result of a conspiracy occurring at the level of the con- 
straints which ensures that the Hamiltonian of the theory is positive definite. When the 
spacetime geometry is spherically symmetric there also exists a quasi-local mass (QLM) 
which is positive and reduces to the ADM mass at infinity in an asymptotically flat geom- 
etry [17]. Attempts to find an analogous quantity when this symmetry is relaxed which is 
also positive have failed. 

In Sect. 4, the QLM of a spherically symmetric geometry is introduced as an integral 
over a spherical surface of fixed proper radius of a spacetime scalar quantity, and, as such, 
an observable of the theory by any reasonable criterion [18]. When an appropriate linear 
combination of the constraints is satisfied the QLM can be expressed as a volume integral 
over the sources. An equivalent expression was derived by Fischler et al. in ref.[17] (see 
also [19]). The QLM thereby provides a very useful first integral of the constraints. 

In Sect. 5 we provide two proofs of the positivity of the QLM when the geometry is 
regular. 

If the dominant energy condition (DEC) ([10]) and the constraints are satisfied posi- 
tivity can be established in a manifestly gauge invariant way [15]. This is achieved remark- 
ably easily by exploiting the optical scalars. We comment on the approach to a singularity 
when the QLM is negative. The second proof is weakly gauge dependent. However, it has 
the peculiar property of permitting us to replace the DEC by the weak energy condition 
and ignore the momentum constraint when we use a linear extrinsic curvature foliation of 



spacetime. 

Both of these positivity proofs arise as simple coroUaries to the existence of appropriate 
bounds on the phase space variables; in the former case an upper bound on the product of 
the optical scalars [15]; in the latter, by the bound on the derivative of the circumferential 
radius of the geometry with respect to its proper radius: — 1 < 9^i? < 1. This bound has 
a simple geometrical interpretation in terms of the embedding of the geometry in flat i?^. 

When the DEC is satisfied and the geometry is regular, additional bounds can be 
placed on the values assumed by the optical scalars, which, in turn, imply a bound on 
the extrinsic curvature. Considering the identification of these variables as the momenta 
conjugate to the spatial metric this is a particularly intriguing result. These bounds assume 
a particularly simple form when the foliation is maximal. It is not clear what role these 
bounds play in the theory. They do not appear to be related directly to the positivity 
of the QLM. It is possible, however, that they will prove to be more fundamental. The 
derivations we provide in Sect. 6 are more economical and the bounds tighter than those 
derived in [15]. 

In Sect 7. we compare the values of the QLM and the material energy. Any reason- 
able defined measure of the gravitational binding energy should always be negative. In 
particular, we demonstrate that the naive definition consisting of the difference between 
the QLM and the material energy is negative when the foliation of spacetime is maximal. 

In Sect. 8 we summarize our understanding of the constraints in terms of the optical 
scalars. This is done by associating with each regular solution of the constraints a closed 
bounded trajectory on the parameter space of the optical scalars. The set of all such 
trajectories can be identified as the phase space of the theory. 

Because of the importance of instantons in the semi-classical approximation, we will 
also occasionally comment on the form of the constraints in Euclidean signature relativity. 

We finish with a summary and an outline of subsequent papers [20] [21]. 

2 THE CONSTRAINTS 

2.1 The Constraints in terms of metric variables 

Initial data for the gravitational field in general relativity consist of a spatial metric 
Qah and an extrinsic curvature tensor Kah which satisfy the constraints [1,2]: 
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K^ -K''^Kab + n = 16np, (2.1a) 

and 

V'Kab - ^aK = -8nJa • (2.16) 

TZ is the three-scalar curvature constructed with gab- The three-scalar p is the material 
energy per unit physical three-volume. The three-vector J" is the corresponding current. 
When the signature of the spacetime metric is made positive definite the sign of the 
quadratic terms in the extrinsic curvature appearing in Eq.(2.1a) is reversed. 

We will examine spherically symmetric spacetime geometries. The only non-trivial 
space-time directions are the radial and time directions orthogonal to the orbits of rotations 
and the geometry can be described by a line element of the form 



ds'^ = -(A^2(r, t) - /^^(r, t))dt^ + 2/3(r, t)dtdr + C{r, tfdr'^ + R^{r, t)dn^ . (2.2) 

The spatial geometry at constant t, we parametrize by two functions C and R of the radial 
coordinate r. A^ and /3 are respectively the lapse and the radial shift. The scalar curvature 
of the spatial geometry is now given by 



—I 



2(rR'\ -R'^-1 



(2.3) 



We introduce the notation ' to denote the derivative with respect to the proper radius i 
defined by d£ = Cdr. When radial derivatives are taken with respect to £, C no longer 
appears explicitly in the constraints. The requirement that this condition be preserved 
under the dynamical evolution of the spatial geometry will determine (3 implicitly. In 
general P will not be zero. 

The other invariant geometrical measure of a spherically symmetric geometry is the 
circumferential radius R. The identification of r with R is the radial coordinate choice 
which is most frequently adopted.* The difficulty, however, is that this identification 



* Indeed, with respect to this choice of radial coordinate and a foliation by Kji = 0, the 
Hamiltonian constraint reduces to an exactly solvable linear first order differential equation 
for C~^. In addition, (3 = 0. 
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breaks down wherever R' = which is the condition that the two-surface of constant r be 
an extremal surface (see appendix A) of the spatial geometry. By comparison, £ increases 
monotonically as we move out from the base of the geometry, insensitive to the formation 
of extremal surfaces (or apparent horizons) so that the identification of i with r is globally 
valid. 

We can write the extrinsic curvature in the form consistent with spherical symmetry 

Kab = riaUbKc + {gab - nanh)Kji , (2.4) 

where Kc and K^ are two spatial scalars and n" is the outward pointing unit normal to 
the two-sphere of fixed r, n" = (£~^,0, 0). With respect to the proper timelike normal 
derivative (A^ = 1 and A^^ = 0), Kab = ilab/'^-, so that Kc = t/ C and Kfi = R/R. Kc is 
also proportional to the acceleration of a radial spacelike geodesic curve on the initial data 
surface. 

The quadratic in Kab appearing in the Hamiltonian constraint can be expressed in 
terms of Kc and K^. The constraints are now given by 

Kr [Kr + 2Kc] - -^ [2 {RR')' - R'^ - 1] = 87rp (2.5) 

and 

K'R + ^{KR-Kc)=4nJ, (2.6) 

where we define the scalar J = J ■ n. All but the radial component of the current three- 
vector J vanish. The only non-vanishing momentum constraint is the projection onto the 
radial direction. 

In any realistic model, matter will be modelled by a field theory, p and J will then 
be cast as functionals of the fields and their momenta. However, for our purposes we 
will suppose that we are given two functions p{£) and J{£) on some compact support, say 
[0, £0] ■ An important fact we will discuss in detail below is that a solution which is both 
asymptotically fiat and non-singular will not exist for every specification of p{£) and J{£)- 
This might happen if an excessive energy, in a sense which will be defined more precisely 
in subsequent papers, is concentrated within a confined region [20,21]. 

We stress that our specification of the sources possesses a spatial diffeomorphism 
invariant meaning. This should be contrasted with the provision of p or J as functions of 
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the flat background coordinate in conformal coordinates, rC = R, with respect to which 
the hne element assumes the conformaUy flat form 

ds'^ = C'^{dr^+r^dn^). (2.7) 

Like the proper radial identiflcation, this system is globally valid. The disadvantage is the 
unphysical nature of the background spatial geometry. Even the simple constant density 
star is not without its subtleties in this gauge despite the fact that the constant density 
is a spatial diffeomorphism invariant. The reason is that the dimensions of the physical 
support of the star is determined in terms of its coordinate dimensions with respect to the 
flat background only after we have solved the constraint. 

In conformal gauge, an appropriate conformal scaling of p is often introduced in order 
to guarantee existence of a solution to the constraints [3]. The result is that one appears to 
be able to sidestep the very singularities we take pains to focus on. While this is flne when 
one is only interested in existence, simply consigning boundary points on the conflguration 
space to inflnity does not help to clarify the physics which underlies the occurence of 
singularities. 

To be fair there is no procedure for solving the initial value problem which is entirely 
satisfactory. Even though the speciflcation of p as a function of i does possess a spatial 
diffeomorphism invariant signiflcance, we have no quantitative notion of the proper volume 
it occupies or, indeed, if such a p can be even consistently specifled until we solve the 
constraints. In the former case, we could, of course, treat V itself as our spatial coordinate. 
This would correspond to the identiflcation 4nR^C = 1. The constraints then provide an 
equation for R (and thus trivially also for C). However, the beneflt we gain is offset by the 
increased non-linearity of the equations. 

2.2 Boundary Conditions 

We are interested in geometries which possess a single asymptotically flat region. It 
is then appropriate to require that the geometry be closed at one end, £ = 0: 

R{0) = (2.8a) . 

In this way, we exclude the possibility that the geometry possess a wormhole to another 
asymptotically flat region or that it does something like degenerate into an inflnite cylinder 
at this end. Local flatness of the metric at this base point also requires that 
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R'{0) = 1 . (2.86) 

A remarkable feature of the constraints is that once we demand that the geometry be 
regular at its base point, this boundary condition is automatically implemented when the 
constraints are satisfied. The only boundary condition we need to impose on R is (2.8a). 
The technical reason for this is the singularity of the Hamiltonian constraint Eq.(2.5) as a 
second order ODE. This is obvious if we rewrite the constraint in the form 

RR" = ^(1 - R'^) + ^Kr[Kr + 2Kc] - 47rR^p. (2.5') 

The right hand side is regular if R is and Kab blows up no faster than R~^. Because R now 
multiples the second derivative the equation must be singular aX R = 0. Once we impose 
the boundary condition (2.8a), however, the requirement that R" also be finite enforces 
Eq.(2.8b) (by convention, we choose the positive sign) and in turn, R"{0) = 0. For a given 
p{£) and J{i), a non-singular solution of the constraint will be unique if it is regular at 

e = o. 

In particular, we will also see that the single boundary condition (2.8a) is sufficient to 
guarantee that spacetime be asymptotically fiat, i? — » £ as £ -^ oo provided the sources are 
distributed on a compact support (or fall off appropriately) and provided the geometry is 
non-singular. 

As an illustration of what might go wrong if i?'(0) 7^ 1, let us compute the three-scalar 
curvature for the spatial geometry described by the line element, (is^ = d£'^+a'^£'^dQ'^ where 
a is some positive constant. In this geometry, -R'(O) = a. If a 7^ 1, the geometry suffers a 
conical singularity at the origin associated with the solid angle deficit, AO = 47r(l — a^). 
This manifests itself in the divergence of 7?. given by 

as the origin is approached. The sign of TZ depends on the sign of a — 1. It is positive 
when a > 1 (a solid angle surplus) and negative when a < 1 (a solid angle deficit). Unlike 
a two-dimensional cone which is fiat away from its apex (7^ = when £ 7^ 0), the conical 
singularity we are considering has a long range field associated with it. In fact, the falloff 
in 7?. is so slow that the space is not even asymptotically fiat. This is the generic behavior 
associated with a conical singularity. Two-dimensional conical structures are exceptional 
in this regard. 
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If both J = and Kab = the hamiltonian constraint gives us that TZ wiU be finite 
when p is. The constraints therefore forbid simple conical singularies (finite R' ^ 1) under 
these conditions.* They do, however, admit more serious cusp singularities (infinite R') 
with a divergence in the traceless component of TZab- If J 7^ 0, however, the constraints do 
not necessarily imply that TZ is finite. This is because a divergence in 7?. can be balanced by 
a divergence in Kab- However, what is true is that R' will always diverge at the singularity 
so that conical singularities cannot occur. The formation of singularities will be discussed 
in sect.5 and in greater detail in II (J = 0)[20] and III(J 7^ 0)[21]. 

2.3 The Constraints in terms of the Optical Scalars 

A remarkable feature of the constraints when the spacetime geometry is spherically 
symmetric is that the constraint equations, (2.5) and (2.6) can be expressed in a symmetri- 
cal form with respect to the optical scalars, defined in terms of the divergence of the future 
pointing and past pointing outward radially directed light rays on the spherical surface of 
fixed proper radius. In appendix A, we show that 

e± = ^{R'±RKR), (2.9a, 6) 

i.e G_(_(G_) is the tangential projection of the sum (difference) of the metric connection 
and the extrinsic curvature tensor. In addition, G_|_ and G_ are canonically conjugate 
variables. In the quantum theory, the G_|_ representation appears to provide a very simple 
characterization of the states which correspond to configurations without apparent horizons 
of the form, *(G+) = 0, if G+ < 0. 

We can invert the defining equations (2.9a) and (2.9b) in favor of the tangentially 
projected two-extrinsic and three-extrinsic curvatures, 

i?'=f(G+ + G_), 

I (2.10a, 6) 

RKn=jie+-e-). 

It is straightforward now to demonstrate that by adding and subtracting an appropri- 
ate linear combination of the constraints, Eqs.(2.5) and (2.6), we obtain two equivalent 
constraints {uj± = RQ±) 



* One could consider a distribution of p which diverges like £ ^ at the origin so that its 
integral over the spatial volume in the neighborhood of £ = is finite. 
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{u+y = - 8nR{p - J) - -^{oo+uj- - 4) + u+Kc 

^^ (2.11a, 6) 

{u-Y = - 8nR{p + J) - 4^(^+^- - 4) - uj-Kc • 

We note that Eq.(2.11b) obtains from Eq.(2.11a) under time reversal, J -^ —J and Kab -^ 
—Kab* In this form, the two constraints are hnear in the 'momentum'. In this sense 
they are the natural 'square roots' of the Hamiltonian. Note, however, that Kc appears 
on the RHS of Eqs.(2.11). The most natural way to treat Kc, in this context, is as an 
independent initial datum specifying some extrinsic time foliation.* We note that these 
equations are simpler than the equations written down by Malec and O Murchadha who 
treat the trace of Kab, instead of Kc, as the independent foliation datum [15]. 

If the geometry is locally flat at the origin, so that Eqs.(2.8a) and (2.8b) hold, and 
the tangential projection of the extrinsic curvature diverges at the origin slower than R~^ 
then 

u+{0) = 2 = u-{0) , {2.12a, b) 

If the geometry is asymptotic flat, in addition, 

lim cu-i- = 2 = lim uj- . (2-13) 

R — >C!0 R — >oo 



* A spacetime argument is presented in appendix B which makes the existence of two 

two such equations more obvious. 

* This is not, however, the usual way to fix such a foliation, which generically will be 

some functional relation of the form, F{Kc, Kji) = 0. 
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3 FOLIATIONS AND SOLUTIONS OF THE MOMENTUM CONSTRAINT 

To fix the foliation of spacetime, we will freeze some homogeneous linear combination 
of the extrinsic curvature scalars. Such a choice is natural because the momentum con- 
straint is itself linear in Kah- Any non-linearity or inhomogeneity in the gauge condition 
would destroy the linear scaling of Kab with J. 

More specifically, let us consider the one-parameter family of foliations, 

Kc + aKR = 0. (3.1) 

We do not consider gauges involving higher spatial derivatives of the extrinsic curvature 
scalars. The foliation (3.1), whenever valid, fixes one (linear combination) of the two 
geometrical momenta at each point. The remaining one is determined completely in terms 
of the intrinsic geometrical and matter variables by solving the momentum constraint. 
This reads 

K'j, + {l + a)^KR = 4nJ, (3.2) 

and is exactly solvable. The solution which is regular at the origin is 

Kr=^I^ dm'+'^J. (3.3) 

Kab vanishes if J = everywhere. In particular, the foliation of Minkowski space by any 
one of these gauges is the standard fiat slicing. If Kr vanishes at any point, then so also 
does Kc so that Kab = there. We note that the boundary condition (2.8a) at the origin 
implies that 



2 + a ' 

in the neighborhood of £ = 0. In particular, Kii{0) = and therefore so also we have that 
Kab = at the origin. 

The slowest acceptable falloff of the extrinsic curvature in an asymptotically fiat ge- 
ometry must be faster than R~^'^ [22]. If a < 0.5 the solution of Eq.(3.2) is inconsistent 
with asymptotic fiatness of the metric in a spatially open geometry. 

If a > 0.5 and is finite, the gauge is globally valid on any regular spatial geometry 
regardless of the support of the initial data, or the presence of extremal or trapped surfaces. 
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It is odd, therefore, that the simplest choice of gauge in the sphericaUy symmetric 
context appears to be the polar gauge * 

Kr = 0, (3.4) 

which corresponds to a — ;► cxo in the parametrization (3.1).t What is most alluring about 
this gauge in a spherically symmetric geometry is that when K^ = 0, the dependence 
on Kab drops out of the Hamiltonian constraint (2.5) which then mimics its form in a 
momentarily static configuration (MSC) Kab = 0. The Hamiltonian constraint can then 
be solved independently of the value assumed by the unfixed extrinsic curvature scalar Kc. 
One particular peculiarity of this gauge is that apparent horizons show up as extremal 
two-surfaces of the intrinsic spatial geometry (see appendix A). Remarkably, in fact, the 
physical content of the model gets encoded completely in this geometry. 

Where then is the snag? To see what price we have to pay, let us examine the solution 
of the momentum constraint (2.6) 

^£ = -4rr^J. (3.5) 

Kc is determined locally in terms of the source. This contrasts dramatically with the 
solution Eq.(3.3) when a is finite where Kc is determined non-locally in terms of J. As 
a tends to infinity, the differential equation (3.2) becomes singular. The support of the 
integrand appearing in the solution collapses in this limit and we recover Eq.(3.5). 

Technically, this is because no derivative of Kc appears in the constraint. The conse- 
quence, however, is that the gauge will break down whenever R' vanishes on the support 
of J. The vanishing of R' signals the development of an extremal two-surface in the spa- 
tial geometry so that the gauge breaks down whenever a current flows across an apparent 
horizon. 

The foliation gauge condition should also fix the lapse. It is easily demonstrated that 
if R' = anywhere on the support of p, the lapse collapses A^(^) -^ in polar gauge 



* What might appear to be the other 'natural' possibility 

Kc = 

corresponds to a = and therefore does not yield a satisfactory falloff. 

t In phase space, this is expressed as He = where He is the momentum canonically 

conjugate to C. 
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— another manifestation of the breakdown of the gauge. Polar gauge clearly does not 
provide a useful description of the physics in the strong field regime (inside matter) we are 
interested in. 

Outside the support of J, Eq.(3.5) implies that Kc = so that Kab = 0. This means 
that the slicing of spacetime defined by polar gauge coincides with the sequence of level 
surfaces of the timelike killing vector in this region. Vacuum spacetime therefore appears 
'static' in polar gauge. This is the optimal exterior form of the metric. In fact, as we 
can see, polar gauge is the unique member of the one-parameter family possessing this 
property. The only obvious shortcoming of finite a-gauges is that they do not provide a 
static slicing outside the support of the matter. 

It is also clear, however, that a static description of spacetime can be approximated 
arbitrarily closely outside the support of matter by letting a be appropriately large. This 
suggests the possibility of constructing a gauge which displays the exterior behavior of 
polar gauge, while at the same time side-stepping its interior shortcomings. What we can 
do is to admit a space dependent parameter a(£) in Eq.(3.1) which tends to infinity outside 
matter. Let a{l) = a + /3(^), such that lim^^oo /3(^) = oo. Then 

Kr = -t:tj—^ dm^+"FJ, (3.3') 

i?i+«F Jo 

where F{(3,£) = exp ( J^ d£(3R'/R). There does not appear to be any gauge (intrinsic or 
extrinsic) which is not tuned artificially by hand which will provide a static description of 
spacetime outside matter. 

The gauge (3.4) is clearly not the only linear combination of the geometric momenta in 
which the Hamiltonian constraint (2.5) mimics the MSC form. The gauge 2Kjc + Kji = 
will also do the job. While this does not appear to suffer from the pathologies of polar 
gauge, it suffers from the shortcoming of producing a slow falloff ~ R~^'^ in Kab outside 
the support of the current. This complicates the asymptotic analysis of the field. As we 
will see, the conventional expression for the ADM mass is no longer valid. 

The existence of two gauges mimicing a MSC is a consequence of the Lorentz signature 
(— , +) of (the metric part of) the supermetric which permits the term quadratic in the 
metric momenta to factorize. These two gauges define the null directions in superspace 
with respect to the supermetric. With respect to a foliation defined by any other linear 
combination of Kc and Kn, extrinsic curvature will show up in Eq.(2.5). 

The two MSC mimicking gauges with a = 0.5 and a = oo define the light cone of the 
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superspace. The admissible gauges constructed using linear combinations of Kc and Kji 
therefore correspond to tangent vectors lying strictly inside the light-cone of the superspace 
metric[16]. The trajectory in the configuration space therefore takes place along spacelike 
directions. This suggests a special role for the light-cone in this mini-superspace. 

The maximal slicing condition K = Kc + 2Kfi = corresponds to a = 2. This is 
the gauge which most readily facilitates the analysis of the constraints in York's conformal 
approach to the full theory and remains the most popular choice among the more formally 
inclined workers in the field. In the spherically symmetric asymptotically fiat context, 
however, this is not such a convincing criterion. Any valid a would appear to offer the same 
reasonable compromise between acceptable asymptotic falloff and non-singular behavior 
in the interior. The remarkable nature of maximal slicing will, however, become evident in 
Sect. 5 within the framework of the optical scalars introduced in Sect. 2. 3 when we specialize 
to initial data satisfying the dominant energy condition. This is not at all obvious in the 
context of the metric variables we have been exploiting in this section. 

We note that when a = 1, the integral appearing in Eq.(3.3) is simply the proper 
volume integral of the radial current scalar J in the interior of £. As we will see in paper 
III, various results simplify dramatically in this gauge.* 

We note that the solution of the momentum constraint requires us to integrate out 
from £ = 0. KjiR^^°' tends asymptotically to the constant value 

)l+a 



diR'+'^J (3.6) 



determined by the current flow. Thus outside the support of the field flow K^ is pro- 
portional to i?~(^+"^. For example, when a = 2, Kab is, up to a constant, the unique 
spherically symmetric transverse-traceless tensor on R^ (see [23]). 

It might appear that we could just as well have chosen to integrate Eq.(3.2) in from 
infinity and to have concluded that outside the support of the field fiow, the solution is 
always K^ = and that the foliation of spacetime is static. The difficulty with this is 
that the resulting solution will be singular at the origin unless the current is fine-tuned 
appropriately. 

Once the momentum constraint has been solved, we substitute (3.1) and (3.3) into 
Eq.(2.5) and solve for R{£) subject to the boundary conditions, (2.8). We defer the details 
to papers II and III. 



* This gauge, like polar gauge, has a particularly simple phase space representation, 
corresponding as it does to the vanishing of the momentum canonically conjugate to R. 
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4 THE QUASI-LOCAL MASS 

An important feature of the constraints when the geometry is sphericaUy symmetric is 
that they possess a first integral which permits the definition of a quasi-local mass (QLM), 
m{t), over a sphere of fixed proper radius which can be expressed as a volume integral over 
the sources contained within that sphere.* To motivate its definition, as well as to make 
a few observations about asymptotic behavior, let us first consider the momentarily static 
data, Kab = 0. We define 

ruK^o = f (l - (R'f) ■ (4-1) 

This should be viewed as a surface integral over the sphere of proper radius £ of a spherically 
symmetric scalar function. When Kab = 0, it is simple to show that the Hamiltonian 
constraint can be cast in the form 

mK^Q=ATxR^R'p, (4.2) 

where m is given by Eq.(4.1) for all values oil. In particular, outside the support of matter 
m'j^^Q = so that m assumes a constant value, TOoo say. This is the ADM mass. If we 
implement regularity at £ = 0, Eq.(2.8), we obtain 

mK=o =4n MR^R'p . (4.3) 

Jo 

Asymptotically, we can now rewrite Eq.(4.1) 

R'^ = l- ^^^ . (4.4) 

R ^ ' 

Thus, as i? ^ oo, i? ~ £ to leading order. The ADM mass in encoded in the next to 
leading order, i? ~ £ — rrioo In^. This is turn permits us to identify a simpler asymptotic 
expression for rrioo'- 



* The 'integr ability' of the system should not be surprising once we identify the me- 
chanical analog of Eq.(2.5) corresponding to the identification of I with time. A generic 
two or higher dimensional model will not be integrable in this way. There has been a fiurry 
of research recently on integrable 'one' and 'two' dimensional models in general relativity 

[24]. 
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moo = lim R{1 - R') . (4.5) 

If Kab does not vanish, the naive generahzation is to replace the quadratic R'^ by the 
square of the spacetime covariant derivative, 

m := ^ (1 - V^RV'R) . (4.6) 

Using the fact that Kji = R/R, this yields the expression 



R^Kl R 

m = 



^ + if (1 - (RT) , (4.7) 



2 2 

which depends, as before, only on initial data. To determine the form of the first integral 
of the constraints analogous to Eq.(4.2), we integrate the the Hamiltonian constraint up 
to £ to obtain 



f'^ 



(l-{R'f)=4n MpR^R'-- MKr[Kr + 2Kc]R^R' ■ (4.8) 

^ ^ Jo 2 Jo 

As before, the LHS coincides asymptotically with rrioo- However, it does not coincide with 

rrioo outside the support of matter. 

Let us add the extrinsic curvature quadratic appearing in the definition (4.7) to the 

LHS of Eq.(4.8). We can exploit the momentum constraint (2.6) to eliminate K^ appearing 

in the integral on the RHS in favor of Kr and J. Thus, modulo the constraints, m satisfies 

m^A-K MR^ [pR' + JRKr] . (4.9) 

Jo 

The RHS is the integral of the scalar, p := pR' + JRKr, over the volume bounded by the 

surface at proper radius £ {dV = 4TrR'^di){see [19]). It is clear that, outside the support 

of matter, to is a constant which we again identify as the ADM mass, nioo- * 

We note that if the extrinsic curvature scalar Kr (as well as Kc ) tends asymptotically 

to zero faster than i?~"^'^, (4.7) reduces asymptotically to the same form (4.5) as Eq.(4.1). 

However, when a = 0.5 the asymptotic form of the surface integral Eq.(4.7) is not (4.1). To 

see this, let us examine the asymptotic 'dependence' of the extrinsic curvature contribution 



* We will see below that the LHS of (4.8), like m, is positive everywhere in any regular 
geometry when spacetime is foliated by an a gauge. In contrast to to, it is even positive 
everywhere in Euclidean relativity. 
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to m on the parameter a. We note that — ^-^ ~ i?^~^" tends to a constant if a = 0.5 and 
diverges if a < 0.5. The latter possibihty was rejected in sect. 3 because it was inconsistent 
with asymptoticaUy flat boundary conditions. We can also see that such a foliation yields 
an asymptotically divergent QLM. In particular, it does not coincide asymptotically with 
moo- 

We noted in sect. 2 that if the geometry is non-singular, and the sources have compact 
support then regularity at the origin is sufficient to force asymptotic flatness. This point 
is clarifled using the flrst integral of the constraint encoded in the deflnition of the QLM. 
From Eq.(4.9) it is clear that if both p and J are compactly supported and R' and Kji 
remain finite, the QLM will be a finite constant outside the support of matter. We also 
noted in the last paragraph that if Kfi falls off fast enough then (4.7) also reproduces 
Eq.(4.4) so that R approaches i in the same way as it does for momentarily static data. 

Only one linear combination of the constraints features in the derivation of (4.9). 
It proves extremely useful to exploit this first integral of the constraints, implementing 
regularity at the origin, in place of one or the other of the constraints. In practice, we 
replace Eq.(2.5) by (4.9) (with m defined by (4.7)). If, in turn, we suppose that spacetime 
is foliated by an a gauge, then we can solve Eq.(2.6) for K^ in terms of J obtaining the 
expression given by Eq.(3.3). 

Note that we have eliminated Kc in going from Eq.(4.8) to Eqs.(4.7)/(4.9) without any 
recourse to a foliation gauge condition. Two properties of the constraints have conspired 
to yield the simple form for |U as a local scalar. The first is that Kc appears linearly in 
Eq.(2.5) and therefore linearly in Eq.(4.8). The second is that it appears undifferentiated 
in Eq.(2.6). In both regards it is unlike K^. There is clearly a conspiracy involving both 
constraints permitting the QLM to be expressed in the simple form Eq.(4.9). 

It is extremely useful to cast the QLM in terms of the optical scalars, G_|_ and 0_. 
We get 



R 

m = - 



(l - \iv+u.) . (4.10) 

The quasilocal mass m is seen to be just the Hawking mass [25]. With respect to these 
variables, Eq.(4.9) assumes the form 

m = 7v diR^[p{uj+ + uj-) + J{uj+-uj-)], (4.11) 

Jo 



or 
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m = TT / MR^ [{p + J)u+ + {p- J)uo-\ . (4.11') 

Jo 



The ADM mass is a spacetime diffeomorphism invariant. In particular, it is independent of 
the fohation of spacetime on which it is constructed. The definition of the quasilocal mass 
either in terms of metric variables, Eq.(4.6) or in terms of the optical scalars shows that 
it is a spacetime scalar though its value does depend on the foliation for finite values of £. 
For each value of £, Eq.(4.9) provides a quasi-local observable of the classical theory. In 
addition, these observables are non-trivial. For whereas rrioo as the effective Hamiltonian 
is trivially conserved, the observables defined by the QLM are not. 

We note that, in general, the spacetime covariant derivative of m can be cast in the 
form [9] 

V^m = ^G'^^e^^ef^^V^'R , (4.12) 

where the notation we use has been defined in appendix B. The Einstein equations can 
now be exploited to recover Eq.(4.9) on projecting (4.12) along the radial direction. The 
evolution of to along the (timelike) normal to the hypersurface, t^, is obtained by projecting 
(4.12) onto t^. Note that the radial pressure will occur on the RHS. 



5 BOUNDS ON THE PHASE SPACE VARIABLES AND THE POSITIVITY 
OF THE QLM 

The most important property of the QLM is its positivity everywhere in any regular 
spherically symmetrical geometry. In this section we will demonstrate how this positivity 
arises as a consequence of bounds on the phase space variables. 

5.1 Positivity of m: The dominant energy condition, and a bound on the product of the 
Optical Scalars 

Let us suppose that the material energy current four-vector is timelike so that it 
satisfies the dominant energy condition (DEC) [10], 



P>^J^a. (5.1) 
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Suppose also that the constraints (2.5) and (2.6), or alternatively (2.11a) and (2.11b), are 
satisfied. Then m is positive everywhere, independently of how we foliate spacetime, if 
the spatial geometry is regular everywhere. Because m coincides with nioo at infinity, this 
provides us with a generalization of the positivity of the ADM mass. 

Because the sources appear explicitly in the manifestly positive combinations, p ± J, 
when the constraints are cast in terms of uj+ and iv-, it suggests that these are the more 
appropriate variables to use when the DEC can be exploited. 

Recall the definition of the QLM in terms of the optical scalars, Eq.(4.10). It is clear 
that the positivity of the QLM is entirely equivalent to the statement 

u;+u;-<4. (5.2) 

This inequality was first derived in [15] but, for completeness, we give a derivation here. 
We note that we can exploit Eqs.(2.11a) and (2.11b) to obtain 



(uj+uj-Y = — Stt-R {u+ + uj-)p — {uj+ — uj-)J — {uj+ + uj-){uj+uj- — 4) . (5-3) 

Note how Kc which appears in both Eqs.(2.11a) and (2.11b) has dropped out of Eq.(5.3). 
This equation is entirely equivalent to Eq.(4.11) with m defined by (4.10). We note first 
of all that the product satisfies the boundary conditions 

uj+uj-{0)=4= lim UJ+UJ-, (2.12c) 

R — >oo 

on account of the boundary conditions, Eqs.(2.12) at the origin and (2.13) at infinity, if the 
geometry is asymptotically fiat. In addition, it must be finite everywhere in any regular 
geometry. And, if the product is finite everywhere, it must possess an interior critical point 
at some finite value of £ if it is not constant. At the critical point, the RHS of Eq.(5.3) 
must vanish. Thus 

LJ+LJ- - 4 = -8nR \p - ^+ ~ ^~ j] , (5.4) 

L ujj^ + uj- J 

unless c<;_|_ + uj- = 0. It is now clear that Eq.(5.2) is satisfied when p > \J\. For if a;_(_ and 
UJ- possess different signs (which includes the case where uj+ + U- = 0) then Eq.(5.2) is 
obviously satisfied. Therefore suppose that they possess the same sign. It is then always 
true that 
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W-i- — UJ- 
CU-i- + UJ- 



<1. (5.5) 



Thus the term appearing in square brackets in Eq.(5.4) is manifestly positive whenever 
Eq.(5.1) is satisfied. This estabhshes Eq.(5.2). 

In those regions of the (a;_|_,a;_) plane where c<;_|_ and U- possess different signs so 
that Uj^U- < 0, not only is m positive but, in addition, m > R/2 regardless of whether 
the constraints are satisfied, or that the energy is positive. In particular, m = R/2 on the 
future and the past apparent horizons. 

We note that the absolute maximum of the product, ujj^uj- obtains at the boundary 
values £ = and £ = oo and it is also the fiat space value. We note, however, the 
corresponding values of m are m{0) = and lim^^oo m = m^ 



"oo- 



5.1.1 Negative QLM and the approach to singularities 

If cu-i-o;- > 4 anywhere, the geometry must possess a singularity. How does this occur? 
If we enter a region in which ujj^uj- > 4 with both u;_|_ and uj- < 0, then when the 
DEC is satisfied, Eq.(5.3) implies that 

{u+u-y>0. (5.6) 

The product U-^-U- monotonically increases. Once m goes negative it decreases montoni- 
cally.* In particular, m cannot recover positive values. The barrier, u-^-iV- = 4, with both 
uj-^- and U- < is therefore semi-permeable. 
In addition, when u+U- > 4, then 

^{u++u.)<-l, (5.7) 

so that R' < —1 and decreasing. Therefore, if the circumferential radius is Rq when 
uj^uj- = 4 we know that the solution must crash, i.e., i? ^ in a finite proper distance 
which is less than Rq from that point. 

What would happen if instead we had a source which did not satisfy the DEC so that 
we entered the region uj+uj- > 4 but with both u)+ and U- > 0? Let us further assume 



* We will see in II and III that it is nonetheless remains finite all the way to any 
singularity. 
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that the source changed its nature so that in this region it did satisfy the DEC. Instead of 
Eq.(5.6), we now have 

{uj+uj-y<0, (5.6') 

which means that the solution is being pushed out of the region uj^uj- > 4 and, instead 
of Eq.(5.7) we get 

^{u;+ + iv-)>l, (5.70 

so that R' > 1. Hence there is no way that the solution can now crash within the region 
uj^uj- > 4. Thus the barrier uj^uj- = 4,a;_|_,a;_ > is also semi-permeable, one can go 
from above down but not up from below so long as the DEC holds. 



5.2 Positivity of m: a-gauges, the weak energy condition and a bound on W^ 

A remarkable feature of foliations of spacetime by the a parameterized gauges is that 
it is possible to (1) relax Eq.(5.1) to the weak energy condition (WEC), viz., p > and, 
(2) omit the momentum constraint yet still establish the positivity of m everywhere. 

The proof is again very simple. However, there is no advantage to be gained by 
exploiting the optical scalars. We return to the definition of m given in terms of the metric 
variables, Eq.(4.7). It is clear that m > whenever (-R')^ — 1- We need therefore only 
show that (-R')^ < 1 under the conditions of the hypothesis and we are done. 

We first note that R' must be bounded in any regular geometry. Because R' = 1 both 
at the origin and at infinity, R' must possess some interior critical point. This will occur 
when R" = 0. The Hamiltonian constraint (2.5) then implies that 

{R'f = 1 - SnpR^ + Kr [Kr + 2Kc] R^ (5.8) 

at such a points. In a MSC, if p > at the critical point then it is certainly always true 
that 

{R'f<l. (5.9) 

Unfortunately, if Kab 7^ this will not generally be true unless Kr [Kr + 2Kc\ < at 
this point. There is, unfortunately, no gauge invariant reason why this should hold. If, 
however, spacetime is foliated by any gauge such that Eq.(3.1) holds at least at the critical 
points of R' , then the third term on the right hand side of Eq.(5.1) is given by 
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KR[KR + 2Kc] = il-2a)Kl, 

which is negative if a > 0.5. This is just the condition defining a globaUy vahd a-gauge. 
This completes the proof that (-R')^ — 1- 

A few comments on the proof: 

We note that the inequahty (5.9) is stronger than the inequahty (5.2). It is clear 
that what the lemma we have proved here is a stronger statement than the positivity of 
the QLM. For it is possible that R'^ > 1 but m > (see fig. (8.1)). In paper III we will 
demonstrate that whereas the converse of the positivity of the QLM is false {viz., m may 
be positive everywhere but the geometry singular) the converse of what we have proved is 
also true: the geometry is regular if and only if — 1 < i?' < 1. 

Because the momentum constraint has not featured in this proof, unlike the first proof 
no control is necessary over the material current such as that implied by the DEC. However, 
in the same way that Eq.(4.11) involves only one linear combination of the constraints, so 
does Eq.(5.3). Thus neither proof requires the full constraints. 

We note that when £ ^ oo we again recover the positivity of the ADM mass. The 
proof is interesting because, unlike the general proof [26], it does not require the dominant 
energy condition to be satisfied. 

Let us suppose that m < somewhere, so that R'^ > 1 at that point. However, if this 
is the case, then Eq.(2.5) (or (2.5')) implies that 

R" <0, 

so that R' is decreasing there. This can only occur if R' < —1. Therefore, if the circumfer- 
ential radius is Rq when m{£o) = 0, then R' < —1. We know then that the solution must 
crash, i.e. i? -^ in a finite proper distance which is less than or equal to Rq from that 
point. 

In fact, we do not even require that m < 0. For, as we remarked before, it is possible 
to have R' < —1 but m > 0. Indeed, it is possible that though R' decreases monotonically, 
m nonetheless remains positive. The catalog of possibilities will be discussed in papers II 
and III. 

We finally recall that a complete specification of the gauge was not necessary in the 
positivity proof provided in this section. However, if we had not specified the gauge 
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everywhere R' < —1 was satisfied, we would not have been able to claim that R' decreased 
monotonically as we proceeded outwards. 

5.3 a gauges and embedding in Euclidean i?^ 

If the inequality R'^ < 1 holds everywhere on the interval [0, £] it is clear that the 
interior geometry can always be embedded as a hypersurface in fiat Euclidean i?^. Thus 
any regular spherically symmetric asymptotically flat three-geometry consistent with the 
Hamiltonian constraint in the gauge Eq.(3.1) can be embedded in i?^. Later we will 
encounter (strongly) singular solutions of the constraints which cannot be thus embedded 
[20,21]. t 

More generally, whenever K°'^Kab — K^ is positive, it is clear from Eq.(2.1a) that the 
scalar curvature, 7^, is also positive when p is. Intuitively, one would expect a spatial ge- 
ometry with a positive TZ to be embedded more readily in a low dimensional flat Euclidean 
space (in the best case, as a hypersurface) than a generic geometry. This is, of course, not 
true of the embedding of such a geometry in Lorentzian flat i?^ which requires that the 
initial data be trivial. 

Analogous statements do not exist in Euclidean relativity where the sign of the ex- 
trinsic curvature quadratic is reversed. Indeed, in the Euclidean theory, p does not even 
possess a definite sign. If p is of the form kinetic energy plus potential energy, the sign 
of the kinetic term will reverse in the Euclidean theory. There is, therefore, no analogous 
positive quasi-local mass result for instantons except when Kab = 0. Indeed, R'^ need not 
be bounded at finite values of i. 

5.3.1 A Universal Bound on i? by £ 

We have seen that the only possible approach to a singularity is through R' < —1. It 
is therefore always true that R' < 1 on any slice defined by an a-gauge. If this inequality 
holds everywhere on the interval [0,£], the proper radius of the geometry will exceed its 
circumferential radius at £. This is simply because then, £ — R = f^^ d£{l — R') > 0. 
Equality only obtains when space is flat. 



^ We note that the lowest dimension into which the Schwarzschild solution can be 
embedded is R^ [27] consistent with the singularity of the geometry at R = 0. 
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6 BOUNDS ON uj+ AND c^_ 

In general, the inequality (5.9) will not be valid, even when the dominant energy 
condition is satisfied, if the gauge is not of the form (3.1). The embedding argument we 
have just considered in Sect. 5. 3 is not gauge invariant. 

In this section we will demonstrate that when Eq.(5.1) is satisfied, both o;^ and o;^ are 
bounded when the full constraints are satisfied. The nature of these bounds is very different 
from that of the upper bound Eq.(5.2) we obtained on the product cu+o;-. The proof that 
such bounds exist proceeds, however, in exactly the same way as the proof of (5.2). Before 
proceeding with the proof, it is useful to recast the constraint equations Eqs.(2.11a) and 
(2.11b) in a form which treats the trace, K, rather than K^ as the independent extrinsic 
curvature scalar which will be fixed by an appropriate gauge condition. Eliminating Kc 
in favor of K we get [15] 






+ 



4 - Au+KR - u+u- 
4 + 4U-KR-U+U- 



(6.1a, 6) 



We will prove that 



l^-i-l < |/t| + ( I 



k\'+A 



1/2 



O 



(6.2a, 6) 



where k = Sup R\K\. We again recall that both uj+ and U- satisfy the boundary conditions 
Eqs.(2.12) at the origin and (2.13) at infinity if the geometry is asymptotically fiat. In 
addition, they must be finite everywhere in any regular geometry. And, if so, they each 
must possess an interior critical point at some finite value of £ if not constant. At the 
critical point, the RHS of Eq.(5.3) must vanish. Thus 



1 r 
4R 



200 



+ 



4 - 4u+K+ 



— c^+o;- - Svr i?(p - J) 



where k+ is the value of RK at the critical point of uj+. Exploiting Eqs.(5.1) and (5.2) we 
have 



u>^ -2u+K+ -4 < 0, , 



(6.3) 



It is now clear that 
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/ X 1/2 / X 1/2 

k+-Mk+P + 4J < u+ < k+ + (\k+\'^ + 4] . (6.4a) 

Similarly 

/ N 1/2 / N 1/2 

-K- - Uk-I'^ + Aj <uj- <~K- + i\K-\^ +4] , (6.46) 

where K- is the value assumed by RK at the critical point of uj-. Now both \k+\ and 
\k-\ are bounded by k which completes the proof of Eqs.(6.2a) and (6.2b). We note that 
we can obtain more transparent (though weaker) inequalities, by further approximating 
(|k|2 + 4)1/2 < 1^1 ^ 2, so that |c^±| < 2 + 2\k\. 

More stringent bounds on a;_(_ and uj- can be extracted from Eqs.(6.4a,b) by using 
k'^ = max RK and k1 = max(— i?i^) to give 

-K*_ - (\Klf + 4] <c^+ <<+ (KI^ + 4) . (6.6a) 

/ \ 1/2 / \ 1/2 

-<-M<|^+4J < UJ- < K*_ + (\k*_\'^ + 4] . (6.66) 

These expressions are useful in the special case where K has a fixed sign. Consider the 
case where K > 0. This gives k'^ — and we get the bounds 

-2 < UJ+ < kI + (\Klf + 4) , (6.7a) 

/ \ 1/2 

-<-(KP + 4) <c^-<2, (6.7a) 



and in the case where i^ < we get 

1/2 



-K*_ - (\k*_\'^ + 4] ' <uj+<2. (6.8a) 

-2<uj-<k*_ + (\k*_\'^ + 4) . (6.86) 

Of course, if spacetime is foliated by an a-gauge, then we can do better still. If J > (<)0, 
then c<;_|_ > (<)a;_. 

We can also place a bound on the sum and difference of uj+ and uj-. We exploit 
Eqs.(6.2a & b) to get, for both the sum and difference 

\R'l\RKR\<^n. (6.9) 
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When K = 0, the former bound coincides with the bound (5.7) in a gauges. If K ^ 0, R'^ 
is still bounded if K is. However, the corresponding spatial geometry will not generally be 
embeddable in R^ . The inequality on the difference has no analog if the dominant energy 
condition is not satisfied. 

It is clear that the bound on the sum can be improved. This is because the bound on 
the product, Eq.(5.2), does not permit a;_|_ and u)- to simultaneously saturate their upper 
and lower bounds. We find 



i«'i^i(" + ^)- 



This is most easily checked using the the graphical representation provided in fig. (8.1). 

The inequalities, Eqs.(6.2a) and (6.2b) come cast naturally in terms of K. Despite 
the fact that the pair of equations, Eqs.(2.11a) and (2.11b), superficially appear simpler 
than (6.1a) and (6.1b) the latter provide the more natural presentation of the constraints. 

A privileged role appears to be played by the maximal slicing of spacetime. Now 
\uj±\ < 2. The bounds Eqs.(6.2a) and (6.2b) now imply the bound (5.2). It might appear 
that when a 7^ 2 the bounds on c<;_|_ and U- are not so useful appearing as they do to 
involve Sup R\K\ explicitly. One can show, however, by bootstrapping on these inequalities 



that that the numerical bound which independent of \K\, \^±\ ^ 2/a/1 — |2 — a], can be 
established in the neighborhood of a = 2 [21]. We now independently possess the bound 
W^ < 1 to estabhsh (5.2). 

7 NEGATIVE BINDING ENERGY 

Once we fix a foliation, a simple measure of the material energy content of the system 
is the (non-conserved) quantity 

M = 47r / MR^ p, (7.1) 

Jo 

also termed the 'bare mass' by ADM [28]. M is a spatial scalar. It is positive and 

monotonically increases with £ when p is positive. It does, however, depend on the foliation. 

It is the QLM which we can think of as the sum of M and a deficit which we tentatively 

identify with the gravitational binding energy Eb associated with the sources in its interior 

m = M + Eb, (7.2) 
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which is independent of the fohation. Its value at infinity is also conserved. 

On the other hand, even though the QLM is positive everywhere when the dominant 
energy condition is satisfied and the geometry is regular, it is easy to see that, in general, 
it will not increase monotonically with I except outside the last apparent horizon. 

To show this, we recall that Eq.(4.10) implies 

m' = 7ri?2 [(p + J)^+ + (p - J)^_] > . (7.3) 

The RHS is clearly positive whenever uj± > (or outside the last apparent horizon) and 
p > \J\ (see Hay ward in Ref.[17]). 

If the initial data possesses an apparent horizon, even though p might be large (so 
that M might also be large if a singularity does not intervene), if p is packed behind 
the apparent horizon nioo can be arbitrarily small. We will examine explicit examples 
in papers II and III. Physically, the material energy is screened by a large gravitational 
binding energy. Because of this the QLM does not provide a very useful measure of the 
material energy. The positivity of the QLM implies that the magnitude of the gravitational 
binding energy can never exceed the material energy in a regular geometry. 

If our understanding is consistent, our definition of the gravitational binding energy as 
the difference between m and M had better be negative at infinity at least. It is surprising 
that, in the maximal slicing, this inequality holds everywhere. M therefore provides a 
global upper bound on m, 

m<M (7.4) 

for all values of £. This is true at a MSC where we can always express the difference [29] 



M-m = 4Tv / diR^p [1 - R'] . 
Jo 



(7.5) 



The right hand side is always positive when p is positive because then R' < 1 everywhere. 
In fact, Bizon, Malec and O Murchadha have shown that we can do much better than this. 
It is possible to place an extremely stringent lower bound on the difference in a MSC [30]. 
We will return to this point in paper II. 
In general, the difference is given by 

M - m = TT / dlR^p [4 - (1 + J/p)uj+ - (1 - J/p)uj-] (7.6) 

Jo 
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= 7r/ MR^p[{l + J/p)2-u+) + {l-J/p){2-u-)] (7.7) 

Jo 

When Eq.(5.1) is satisfied and when iiT = 0, so that |a; ± | < 2, the difference is positive. 
In general, with i^ 7^ it is difficult to see how the positivity of M — m will hold. This is 
an open question worth settling. 
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8 THE {u+,u-) PLANE 

We can exploit the {u+,uj-) plane to represent solutions to the constraints of spher- 
ically symmetric general relativity. We cast the constraints in the form (2.11a & b). To 
solve these equations we must supplement them with the equation which anchors a;_|_ and 
uj- to R, 

R' = {u++U-)/4:. (2.11c) 

On the right hand side of (2.11a & b) appear three additional functions; the material 
sources, p{£) and J{£) and the extrinsic curvature scalar, K £,{()■ 

Our approach has been to specify p{l) and J(£) on some compact interval [0, Iq] 
consistent with the dominant energy condition, (5.1), though we have seen that it is possible 
to relax this condition to the weak energy condition under certain circumstances. It is we 
who decide whether they do or don't satisfy the energy condition — this is not something 
we derive. 

We could specify K^ as some function of L However, intuitively, extrinsic curvature 
should respond to to the flow of matter J. In particular, if we were to do this, in the 
absence of sources, we would flnd ourselves foliating Minkowski space non-trivially. It is 
therefore not really appropriate to treat K^ the same way as p or J. What we do is 
specify some foliation gauge appropriate to the topology under consideration, (3.1) say. 
This permits us to eliminate Kc in favor of Kn = {uj+ — uj-)/AR. The right hand sides of 
Eqs.(2.11a &b) now involve only the functions R, a;_|_ and U- and the two functions p{t) 
and J{t). 

The resulting three coupled ordinary flrst order differential equations can now be 
solved subject to boundary conditions on 0;+, U- and R appropriate to the topology, in 
our case these boundary conditions are (2.12a & b) and (2.8a). 

Each solution of these equations will deflne a trajectory F* = (i?(£), a; _(_(£), a;_(£)) on 
the space of triplets, (i?, a;+, a;_). What is remarkable is that no essential information is 
lost by limiting ourselves to the projections of these trajectories onto the {us^.uj-) plane. 

Whenever the dominant energy condition holds all regular asymptotically flat trajec- 
tories are bounded on the (a;+, a;_) plane by Eqs.(6.2a) and (6.2b) as well as the positivity 
of the QLM, Eq.(5.2). Because these inequalities are independent of R when cast with 
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respect to these variables, they can be represented in the projection. The region E of the 
plane in which these inequalities are simultaneously satisfied is indicated in fig. (8.1). Reg- 
ular, asymptotically fiat solutions are confined to E. Any trajectory which strays outside 
E is necessarily singular. Not only can it not reenter E, the trajectory must run off to 
infinity on the (aj+,a;_) plane at some finite value of I. The details will be discussed in 
papers II and III. 

If the gauge is maximal, this is a square region with vertices (2, 2), (2, —2), (--2, —2) 
and (—2,2). The inequality Eq.(5.2) is a consequence of the other two inequalities. In 
general, however, (5.2) bites out two discs from the square defined by the other two in- 
equalities. In this paper, the bounds we derive on Uj^ and U- characterized by the number 
O do depend on the trajectory itself unless K = Q. derived — not something we put in 
by hand. However, the fact that this bound changes from trajectory to trajectory when 
iiT 7^ is not entirely satisfactory. We can do better. In paper III we show that in any of 
the gauges (3.1), O can in turn be bounded by a universal numerical constant that depends 
only on the parameter a appearing in the gauge. 

The boundary conditions, (2.12) and (2.13) imply that each non-singular trajectory 
must both begin and end on the point, (2, 2). Thus, physical non-singular initial data can 
be identified with bounded closed curves, Fq, in E each of which contains the point (2, 2). 
In general, trajectories can intersect themselves any number of times. The configuration 
space of spherically symmetric general relativity can be identified as the space of all such 
bounded closed trajectories. We note 

1. Vacuum, fiat data corresponds to the zero trajectory, V = (2, 2) for all £. 

2. Initial data which do not possess apparent horizons correspond to trajectories which 
lie in the upper right hand quadrant, uj+,uj- > 0. 

3. Spatial geometries which do not possess extremal surfaces, R' > 0, correspond to 
trajectories which lie to the right of the principal (negatively sloped) diagonal, R' = 

4. Moment of time symmetry initial data correpond to trajectories lying on the positively 
sloped diagonal, Kfi = 0. The extremal surface condition coincides with the apparent 
horizon condition. 

If spacetime is foliated by an a gauge, in addition, we have the inequality, (5.7), 
satisfied by R'. If a 7^ 2, this reduces the range of the allowed trajectories further, reducing 
E to the hexagonal E^, illustrated in fig. (9.1). 

5. In an a gauge, when J > (J < 0) everywhere so also is K^. Thus the trajectory 
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lies below (above) the diagonal, Kji = 0. 
6. The existence of an extremal surface does not necessarily imply that of a future ap- 
parent horizon. Why this is so is clear when J > in an ct - gauge. Conversely, the 
existence of a future apparent horizon does not necessarily imply the existence of an 
extremal surface. 

A small loop in the neighborhood of the point (2, 2) corresponds to almost flat initial 
data. The length of such a trajectory will typically be close to zero — the flat data value. 
A larger loop, crossing a;+ = has a length bounded from below by 4. It would appear 
that the length of a trajectory corresponds somehow to how far the initial data is from 
vacuum flat initial data. This criterion does not, however, require the trajectory to venture 
far from the point (2,2). For example, a trajectory might wiggle about so much that it 
possess an arbitrarily large arc length even though the distance of maximum excursion 
from (2, 2) is never large.* In paper II, we flnd better criteria to characterize the distance 
of the initial data from flat space by deflning a norm on the space of initial data which 
consigns singular initial data to inflnity — something the naive idea proposed here fails to 
do. 

9 CONCLUSIONS 

This article has been devoted to an examination of the constraints in spherically 
symmetric general relativity with the goal of identifying the physical degrees of freedom. 
With this goal in mind, we found that it was useful to exploit not only the traditional 
canonical description of the phase space provided by the metric variables, {gab, Kab), but 
also the optical scalar variables, (a;+,a;_). An intriguing feature of the optical scalars is 
the possibility of casting the constraints in the linear form (2.11a and b). Working with 
the appropriate set of variables, we could focus in on different properties of the phase 
space; gab to describe spatial metric properties; Kab to describe the foliation; us^ and U- 
to describe the light cone. Solutions are represented as trajectories on the (c<;_|_,a;_) (or 
{R',RKr)) plane as described in Sect. 8. 

Our approach distinguishes between properties which are spacetime diffeomorphism 



* We note that the points extremizing the length of the excursion from flat space, 
{ujj^ — 2)^ + {uj- — 2)^, are the natural analogs of the points at which R" = at a moment 
of time symmetry. 

36 



invariant — independent of the foliation, and those which are not. 

When we do fix the fohation, we take particular care to introduce gauges which are 
global and which do not break down at either minimal surfaces or at apparent horizons. 
We introduced a one-parameter family of linear extrinsic time foliations of spacetime in 
Sect. 3 which includes both the polar gauge and the maximal slicing condition. It turns out 
that only a subset of these give reasonable asymptotic falloffs to the initial data and these 
gauges are those bounded by the null directions of the superspace metric. These gauges 
are the natural asymptotically fiat foliations of spacetime. If the spherically symmetric 
model is any indication, the null directions in superspace may provide a guide towards 
identifying the natural gauges in less trivial models. 

In spherically symmetric general relativity it is easy to identify a spacetime diffeomor- 
phism invariant quasi-local mass, which coincides with the Hawking mass, and to write 
an expression which relates this quasi-local mass to an integral over the sources when 
the constraints are satisfied. The remarkable feature of the QLM is its positivity when 
the geometry is regular and the material sources satisfy certain very reasonable energy 
conditions. The positivity of the QLM appears to suggest the negativity of a physically 
realistically defined binding energy. 

We found that we could prove the positivity of the QLM under two different sets of 
assumptions on the initial data. One of the proofs is gauge invariant, the other is not — 
relying, in addition, on the implementation of a valid a. - gauge. However, whereas the 
former requires that matter satisfies the DEC (as one would expect), the latter does not 
— requiring only that matter satisfy the WEC. This is not to say that the conjunction of 
the WEC with an a - gauge is equivalent to the DEC. 

Underpinning the positivity in either case, are various bounds on the canonical vari- 
ables. When the former (latter) set of conditions mentioned in the preceeding paragraph 
are satisfied, cu+cu- < 4 (i? ^ < 1) everywhere in any regular solution to the constraints. 
This bound on B! has a simple geometrical interpretation in that it demonstrates that each 
of the spacelike hypersurfaces we obtain as solutions to the constraints can be embedded 
in fiat i?^. In fact, given the spherical symmetry, they can be described simply by curves 
in i?^. Of course, this visualization ignores the nontrivial extrinsic curvature. 

There is no analogous bound on K^ if the DEC is not satisfied. One might be forgiven 
for overlook even attempting to search for such a bound because typically one does not 
expect momenta to be bounded. What is remarkable is that when the DEC is satisfied, K^^ 
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is bounded. This is not at all obvious using the metric variables. The way one proceeds 
is to establish that both uP\^ and uPl are bounded when the DEC is satisfied. It is then a 
simple corollary that K^ is also bounded. Clearly the spacetime light cone structure on 
the initial hypersurface, though heavily disguised in the metric description of the initial 
data, is encoded in the constraints. While the guiding principle behind the discovery 
of the bounds on uj^ and U- may be the positivity of the QLM, these bounds are of a 
fundamentally different nature to that on the product, o^+a;-, which features in the proof 
of the positivity of the QLM. We note that the form of these bounds is also very different 
from that of the bounds separating geometries with apparent horizons from ones which 
do not. The existence of bounds on the canonical variables (or their gradients) appears 
to be a very fundamental feature of the spherically symmetric theory, undermining one's 
confidence in any naive quantization of this model which does not take them into account. 
It is likely that analogous bounds will exist in the full theory. Their discovery is a challenge 
for the future. 

What we have learnt about the configuration space of spherically symmetric general 
relativity will be built upon in future papers. In particular, we will show that two ingredi- 
ents, the description of solutions as trajectories on the (cu-i-, a;_) plane and the quasi- local 
mass provide extremely useful practical tools, not just abstract constructions. We can 
exploit the QLM to characterize the behavior of the spatial metric and the extrinsic curva- 
ture in the neighborhood of generic singularities. The only singularities that can occur in 
a spherically symmetric geometry do so because K returns to zero. Generically this will be 
accompanied by a divergence of both B! (towards minus infinity) and of K^. No singular 
geometries terminate in the region, E, indicated in fig. (8.1). Where they terminate will 
depend in an essential way on the value of the QLM in the neighborhood of the singularity. 

A very brief outline of subsequent papers is: 

In II, we will examine the solution of the Hamiltonian constraint at a moment of time 
symmetry. In this simplified setting we can gain useful clues as to how best to characterize 
the configuration space of the theory. In III, we extend this analysis to J 7^ 0. 
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Appendix A: Optical Scalars, Apparent Horizons and Extremal Surfaces 

Let us consider a closed two-dimensional spacelike surface, S. Each point on S pos- 
sesses two mutually orthogonal spacetime unit vectors which are normal to S. One of 
these, t^ say, may be taken to be timelike and future directed. The other vector, n^ say, is 
then spacelike and we choose it to be pointing outwards. This choice is clearly well defined 
up to a Lorentz boost in the normal tangent space. Alternatively, we can always choose 
two null normal vectors, one outward directed and future pointing, k^ say, the other also 
outward directed but pointing into the past, k'^ say. 

We can represent these null vectors 

k'^ = ±t^ + n^. (a.l) 

/c^ is obtained from /c_|_ by reversing the sign of t^. The divergence of /c^ on S is now given 
by 

e± = {g'"' + tH- - n^n-)W^{±K + n,) ■ (a.2) 

Let us suppose that S is embedded in some spacelike hypersurface S with normal vector 
t^. Now the normal vector to 5" in iS is also clearly normal to t^ in spacetime. With 
respect to Gaussian normal coordinates for spacetime adapted to 5, t^ — (1,0,0,0) and 
vP = Q = riQ. Thus we can express 

©±=(^"'-^"^')V„(±4 + ri6) 

(a.3) 

where gab and Kab = ^atb are respectively the spatial metric and the extrinsic curvature 
of E. The second term is the trace of the projection of Kab orthogonal to 5". In particular, 
we observe that G+ and G_ are completely described by the initial data, (gab^ Kab) on 
the spacelike hypersurface in which we have embedded 5". We note that a change in the 
prescription of E will change G+ by a boost factor, 7 say, and G_ by the factor 7"^. Thus 
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neither G_(_ nor G_(_ is a spacetime scalar. However, their product G_G_|_ is. This is the 
product which occurs in the quasilocal mass formula introduced in Sect. 5. Clearly, m is 
also a spacetime scalar. 

A future (past) trapped surface is a closed two-dimensional spacelike surface on which 
the divergence of future (past) outward directed null rays is negative. A future (past) 
apparent horizon is the outer boundary of such trapped surfaces. The appearance of 
a future (past) apparent horizon signals (Penrose Theorem[10]) gravitational collapse to 
form a black hole (initial conditions which could only have evolved out of a state which 
possesses a singularity). * 

The condition 

V ■ n = (a.4) 

is the condition that the closed two-dimensional spacelike surface S be an extremal hy- 
persurface of S. The apparent horizon coincides with an extremal (minimal) surface in a 
MSC, regardless of the material content of the theory. 

In the spherically symmetric model, the normal spacelike vector to the spherically 
symmetric two-dimensional surface of circumferential radius R is given by n" = ^(1, 0, 0) 
and 

Van" = 4^'. 



In addition, 



^gab _ n-n^)Kab = 2Ki 



This term therefore not only vanishes in a MSC but also when K^ — 0. In general, 
however, it will not. We can write 

G± = |(i?'±i?Kfl). (a.5) 

Recall that with respect to the proper timelike normal, RKji = R. We can therefore 
alternatively write 

-Tt 

In flat space foliated by flat spacelike hyper surfaces, R — £ independent of t, and Kab = 0. 
Thus 

e±R = 2, 



We note that there exists no analog of an apparent horizon in a Euclidean geometry. 
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for all R. 

A minimal surface in the spatial geometry does not necessarily correspond to any 
physically significant locus of points on the spatial geometry. However, it does imply that 
the geometry possesses either a future or a past apparent horizon. Even if an apparent 
horizon is not present on the initial spacelike surface, as the system evolves an apparent 
horizon might form. One of the nice things about the identification of the radial coordinate 
with the proper radius is that it is insensitive to the formation of minimal surfaces or 
trapped surfaces. This should be contrasted with Schwarzschild coordinates, which even 
if globally valid on the initial surface, will not necessarily remain so. 

Appendix B: Spacetime Approach to Eqs.(2.11a &; b) 

Another route to the derivation of Eqs.(2.11a) and (2.11b) which makes them, perhaps, 
more obvious is the following spacetime approach suggested to the authors by the referee. 
We note that the constraint equations (2.1a) and (2.1b) are equivalent to the projected 
Einstein equations {G'^'^ is the Einstein tensor) 



G^,r = 87rT^,r , (6.1) 

where t^ denotes the future directed unit normal to the hypersurface, S. We exploit the 
'radial' Einstein equation [9] 

h'^^h^.V^VpR = ^ V, - 47ri?T«^e«^e^, , (6.2) 

where h^^ is the r — t part of the spacetime metric, m is the QLM 

discussed in sect. 4, and e^^, is a 2-form associated with the surfaces orthogonal to the 
orbits of the rotation group. If n^ be the radial tangent to 5, then e^^, = 2t^^n^,y By 
contracting (b.2) with n^ we obtain a set of equations equivalent to (b.l). The Hamiltonian 
and momentum constraints obtain by projecting the resulting equations onto t" and n" 
respectively. If, however, we project onto the two linear combinations /c^ defined by 
Eq.(a.l), we get 
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kln^Vc^VpR =^T AnRT^pklt^ , (6.3) 

a set of equations equivalent to Eqs.(2.11a) and (2.11b). 

Figure Caption 

fig.(8.1) The (a;+,a;_) plane. Regular asymptotically fiat solutions are confined to the 
region, E, bounded by the closed union of line and arc segments AB,BC,CD,DE,EF, and 
FA. When Eq.(3.1) is satisfied, this is reduced to the hexagonal region E^, bounded by 
the closed union of hue segments, A'B',B'C,CD' , D'E',E'F, and FA'. Both a regular 
trajectory Fq, and a singular one F are illustrated. 
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